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Abstract 

Linear alkylbenzene (LAB) is adopted to be the organic solvent for the 
Jiangmen Underground Neutrino Observatory (JUNO) liquid scintillator de¬ 
tectors due to the ultra-transparency. However the current Rayleigh scatter¬ 
ing length calculation disagrees with the measurement. The present paper for 
the hrst time reports the Rayleigh scattering of LAB being anisotropic and 
the depolarization ratio being 0.31 ± 0.01(stat.) ± 0.01(sys.). We proposed 
an indirectly method for Rayleigh scattering measurement with Einstein- 
Smoluchowski-Cabannes formula, and the Rayleigh scattering length of LAB 
is determined to be 28.2 ± 1.0 m at 430 nm. 
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1. Introduction 

The Jiangmen Underground Neutrino Observatory (JUNO) is a multi¬ 
purpose neutrino experiment designed to determine neutrino mass hierarchy 
and precisely measure oscillation parameters by the medium baseline vacuum 
oscillations of the reactor antineutrinos oca- It’s located in the 700 m deep 
underground laboratory at Jiangmen, China. 

The design of JUNO is based on the principle with highly purihed liquid 
scintillator (LS) as the central detector surrounded by ten-thousands of pho- 
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tomultiplier tubes (PMTs) and tons of ultra-pure water outside as an external 
shield. The central detector is a 20,000 tons of linear alkylbenzene (LAB) 
based LS in a spherical vessel with a diameter of 34.5 m, and the energy reso¬ 
lution is designed to be 3%/y/ii^(MeV) corresponding to 1, 200 photoelectron 
(p.e.) collected by PMTs per MeV. The antineutrino is detected via the in¬ 
verse beta decay reaction Vf>+p^e'^ + n. LAB scintillation light generated 
from positron annihilation and neutron capture is red-shifted by the primary 
and secondary wavelength shifter [3], 2,5-diphenyloxazole (PPO) and pbis[2- 
methylstyryljbenzene (bis-MSB), and has to transverse through the whole LS 
vessel before arriving at the photomultipliers (PMTs). As a consequence, the 
transparency of LAB to the shifted scintillation light (430 nm for instance) 
is essential [1] . This has brought a great challenge to both of technical LAB 
purihcation and optical parameter measurement. 

The event location and energy reconstruction are determined by the op¬ 
tical modeling which includes Rayleigh scattering length, absorption length 
and attention length [5]. However, the absorption length which describes 
the energy of photon being absorbed to heat is difficult to measure directly. 
An indirectly method calculating from the attenuation length and scattering 
length is proposed[6]. Currently the attenuation length of LAB at 430 nm 
has been measured[7l IH [9], while the scattering length has a difference from 
40 m HDl to calculated 30 m m- These motivate a precise measurement of 
LAB Rayleigh scattering length at the scintillation wavelength of LS. 


2. Rayleigh scattering 

Rayleigh scattering, developed by Rayleigh in 1899 |T2], describes the 
light elastically scattering off the medium molecules. For gas state, this the¬ 
ory was successfully applied to independently isotropic molecule, and mod¬ 
ified by Cabannes with introducing a depolarization ratio to describe the 
anisotropy of molecules. For liquid state, due to the strong interaction ef¬ 
fects between molecules, Einstein and Smoluchowski proposed scattering to 
be caused by the random motion of molecules which leads the fluctuations 
of density and the dielectric constant. The Rayleigh length of liquids can be 
described by Einstein-Smoluchowski-Cabannes (ESC) formula [TT] : 
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Here A is the wavelength of scattered light, n is the refractive index, k is the 
Boltzmann constant, T the absolute temperature, /?t the isothermal com¬ 
pressibility and 5 is the depolarization ratio. For JUNO experiment, the 
temperature of the LS detector will be controlled at T = 20 ± 1°C and the 
quantum efficiency of PMTs are optimized at A = 430 nm. Recently, the 
Pt of LAB at three temperatures over 4 to 23°C has been measured by the 
vibrating tube method na. The refractive index n of LAB from the same 
batch in the range between 400 nm and 630 nm has been reported m- 

The depolarization ratio 6 can be measured at 6 = 90° scattering angle 
with a vertically polarized incident beam according to [mui]: 
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The subscripts designate the components analyzed in the scattered beam, v 
and h are vertically and horizontally polarized portion, respectively. If we 
dehne depolarization ratio fraction / = Ih/h, we can get S = 2 //(/ -|- 1 ). 


3. Experimental setup 

This experiment is designed to measure the depolarization ratio ^ 90 ° . 
The schematic setup is shown in Fig. The light source is a Pico-Quant 
LDH pulsed laser diode with adjustable output up to 10 mW at 40 MHz rep¬ 
etition rate. Its wavelength is 405 nm instead of 430 nm because it’s the only 
available pulsed laser we have which is the closest one to 430 nm. The beam 
from the laser has a divergence of 0.32 mRad, producing a beam diameter at 
the sample cell less than 1 mm. The incident beam is collimated with two 
variable apertures, as well as vertically polarized (polarization< 1 %) with a 
Glan-Thompson polarizer. Fluctuations in the intensity of the light source 
are monitored by use of a reference PMT-I (Hamamatsu R2083) with a grey 
hlter in front to block a portion of the incident light. The samples are held 
in a quartz cuvette with 5-cm path lengths. The scattered light from the 
sample is then polarized by a Glan-Laser calcite polarizer with 100, 000 : 1 
extinction ratio which is installed on a motorized rotation mount with 1 " pre¬ 
cision. For this measurement the selected polarization component is either 
vertically or horizontally to the plate dehned by incident beam and scattered 
light. A Hamamatsu R1828 PMT-H with single photon response capability 
is installed after to count the number of scattered photons. The performance 
of PMT-H is monitored with a blue LED. Both laser diode and LED are 
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trigged by a pulse generator. This experiment is held in a dark room with 
room temperature controlled at 23 °C. 


laser 



Figure 1: Sketch of the experimental set-up 


The signals from both PMTs are recorded by a CAEN DT5720 waveform 
digitizer with 4nm time resolution (250 MS/s) and a total range of 2V at 
a 12 bit resolution. The trigger of data acquisition (DAQ) system is shared 
from the same pulse generator, and it’s set to be IkHz with a 4/is time 
window. The trigger for laser diode is set to a burst mode, 60 pulses with 
5 ns time width are sent in this time window. Both of the signal from PMTs 
(usually PMT signals are less than 20 ns) and the thermal noise or radioactive 
background are recorded by DAQ. 

4. Photoelectron counting 

The intensity of scattered light is estimated by counting the number of 
p. e., then the depolarization ratio fraction / can be written as: 


OO 
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However, the horizontal and vertical portion of scattered light are measured 
separately, which requires a careful consideration for normalizing the experi¬ 
mental conditions such as the stability of laser diode, DAQ data taking time 
period, etc. 
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We propose a new method for estimating the fraction /. Assuming a 
Poisson distribution for the photon scattered and the p.e. number leaving 
the photocathode, and taking the number of incident photons to be N, we 
can rewrite this to: 
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Here /x is the expected value for Poisson distribution. 

Due to the long Rayleigh scattering length, the expected value /i would 
be rather small. The output of the PMT could be altered by the dark noise 
which causes a number of random coincidences are detected. Therefore, in 
order to have the random coincidences contribution at the level of 1% it’s 
necessary to have /i > fdark ■ T'gate/O.Ol. Here fdark is the frequency of dark 
noise which is 4K for PMT-H, and Tgate is the ADC gate which is 40 ns. It 
gives jj, > 0.016 to have a negligible contribution of the dark noise spectrum. 

For a smaller n < 0.016, the observed Poisson distribution includes two 
parts. First one is the scattered photon response contribution, and the second 
one is the dark noise of PMT coming from the thermionic emission from 
the photocathode or radioactive background. The latter can be assumed to 
be a Poisson distribution as well. These two Poisson distributions are not 
correlated. Thus the fraction / is: 
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The superscripts mean the observed and the dark noise background Pois¬ 
son distribution. Since the PMT dark noise spectrum is not related to the 
polarization angle, the and should be at the same level. 

To estimate the expected value /i, the PMT charge spectrum is fitted 
with the following function: 
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Here Nq is a normalization factor, P(n; fi) is the probability of Poisson dis¬ 
tribution with mean value jj, for different contribution of n photons. The 
description of the electronics noise function fnoiseix), single photon response 


5 



function fsep{x) and multiple photon response function can be found 

somewhere at na. 

The charge response of the PMT-II for a low intensity light has been 
studied. The different fi of p.e. is achieved by tuning the applied voltage 
on the LED shown in Figj^ The charge spectrum taken with two different 
mean p.e. number are given in Figj^ On the left plot the single photon 
response including an exponential part and a Gaussian one is shown. Here 
the exponential part is to describe the thermionic emission contribution with 
a smaller energy or under unfavorable angles of incidence. The Gaussian part 
is the single p.e. response with the mean value Xi and the standard deviation 
(Ti. The contribution of 2 and 3 p.e. can also be seen. On the right plot the 
multiple-photon response is given. 

5. Result 

The results reported here include two samples. The hrst one is LAB, and 
the PMT spectrum is shown in Figj^ The dots are observed data spectrum 
and the red line is htted spectrum. The shaded histogram is the dark noise 
spectrum and the dashed line is the ht result. The measured depolarization 
ratio fraction is ftAB = 0.186±0.006(stat.), corresponding to Slab = 0.314± 
0.014(stat.). 

The second one is benzene with 99.5% purity. As one of the most com¬ 
mon organic solvents, benzene is deleterious and highly volatile. Gareful 
consideration has been done during this measurement. The depolarization 
ratio fraction is measured to be /sen = 0.287 ± 0.004(stat.), which gives 
Sbeti = 0.44 ± 0.01(stat.). This is agreed within 3 a with the measurement 
made by Mossoulier which is about 0.42 |16j . 

6. Systematic error 

The systematic uncertainty of this measurement has been checked from 
the following respects: 

The wavelength of incident laser beam has been measured with Ocean 
Optics 4400 spectrometer, and Gaussian htted results show the wavelength 
of laser is 405 ± 0.66 nm. The systematic error from this contribution is 
negligible. 

The stability of laser is monitored with PMT-I. The intensity huctuation 
of the incident laser beam is about 0.1%. As discussed in section the 
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expected value /i is not related to the intensity of the incident laser, this 
fluctuation will not be taken into account for the systematic error estimation. 

During the measurement procedure, the cuvette will be washed and air 
dried for couple times. Even though the light path will be adjusted again 
after cleansing, the position of the cuvette will still be slightly different. The 
systematic from this difference is considered by measuring the absorption 
spectrum of cuvette with a ultra-violet spectrophotometer (UV-Vis). By 
changing the position of cuvette in the UV-Vis, the difference of the ab¬ 
sorption value at 405 nm is about 1.7% and this is taken as the systematic 
error. 

The acceptance correction for PMT-II is estimated by introducing a slit 
function. This slit function S{6) can be approximately written as a function 
of the inner open-angle a and the outer open-angle /3 HZ]. In our measure¬ 
ment a = 0.0598°, and (3 = 0.145°, this correction is estimated by calculating 
the integral S{6) ■ [5 — <^ 9 o ]/<^90 ■ d{9), and it shows a negligible contri¬ 

bution from this correction. 
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The precision of polarization angle is guaranteed by a high precision rota¬ 
tion mount with ±0.1° resolution. For the polarization of the incident beam, 
this systematic error is estimated by rotating the polarizer from —1° to 1°, 
and the difference is about 2.2%. For the polarization of the scattered light, 
the systematic error is estimated by taking this uncertainty to the depolar¬ 
ization ratio calculation, and the largest offset is about 1.0% which is taken 
as the systematic error from this contribution. 

The total systematic error is about 2.8%, and the details are given in 
TableHl 


7. Discussion and conclusion 

The depolarization ratio Slab at 405 nm is measured to be Slab = 0.31 ± 
0 .01(stat.) ± 0.01(sys.), here the hrst one is statistical error and the sec¬ 
ond one is systematic error. According to the benzene measurement |16j . 
SBen has a negligible difference between 405 nm and 430 nm. Assuming Slab 
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Table 1: Systematic error 


Laser wavelength 

negligible 

Polarization of incident beam 

2 . 2 % 

Uniformity of cuvette 

1.7% 

Acceptance 

negligible 

Polarization of scattered light 

1 . 0 % 

Total 

3.0% 


has the same value at 430 nm and taking the refractivity[TT] and isothermal 
compressibility [13], the Rayleigh scattering length of LAB is about 28.2 ± 
1.0 m at 430 nm at room temperature. Here the LAB sample used in these 
measurements are purihed from a set of steps including distillation, AI 2 O 3 
sorption, etc. The same purihcation technique will be used for the JUNO 
LAB mass production. 

The attenuation length of the same LAB sample at 430 nm has been mea¬ 
sured to be about 19.44±0.62m[9|. Thus the absorption length is 62.5±7.9m, 
here the error is calculated from error propagation. The simulation shows 
that the energy resolution can reach the design goal with the parameters 
given above [IB]. 

Rayleigh scattering length, as one of the important optical parameters, 
describes the probability of photons scattered after transmitting a distance 
in medium. The shorter Rayleigh scattering length is, the higher probability 
photon being scattered. This increases the complexity of event vertex and 
energy reconstruction. To simplify this problem, it’s better to have a long 
Rayleigh scattering. Our measurement of Rayleigh scattering length is closer 
to the calculated value reported in m than the direct measurement in [TO] . 
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channels 



Figure 2: The charge spectrum taken with two different mean photoelectron number. The 
dots are data, the red line is the fit result, and dashed lines are the contributions of each 
component. 




Figure 3: The charge spectrum for vertical (left) and horizontal (right) portion of scattered 
light. The dots are data and red line is fitted spectrum. The shaded histogram is the dark 
noise spectrum and the dashed line is the fit result. 


11 








































